Draft version February 5, 2009 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



OBSERVATIONAL SIGNATURES OF SIMULATED RECONNECTION EVENTS IN THE SOLAR 

CHROMOSPHERE AND TRANSITION REGION 

L. Heggland 

Institute of Theoretical Astrophysics, University of Oslo, P.O. Box 1029, Blindern, N-0315 Oslo, Norway 



o 
o 

(N 
X) 

in 



6 



> 

o 

O 

o 



X 



B. De Pontieu 

Lockheed Martin Solar and Astrophysics Laboratory, 3251 Hanover Street, Org. ADBS, Building 252, Palo Alto, CA 94304 

V. H. Hansteen^ 

Institute of Theoretical Astrophysics, University of Oslo, P.O. Box 1029, Blindern, N-0315 Oslo, Norway 

Draft version February 5, 2009 

ABSTRACT 

We present the results of numerical simulations of wave-induced magnetic reconnection in a model 
of the solar atmosphere. In the magnetic field geometry we study in this article, the waves, driven by 
a monochromatic piston and a driver taken from Hinode observations, induce periodic reconnection 
of the magnetic field, and this reconnection appears to help drive long-period chromospheric jets. By 
synthesizing observations for a variety of wavelengths that are sensitive to a wide range of tempera- 
tures, we shed light on the often confusing relationship between the plethora of jet-like phenomena 
in the solar atmosphere, e.g., explosive events, spicules, blinkers, and other phenomena thought to be 
caused by reconnection. 

Subject headings: magnetic fields — MHD — Sun: chromosphere — Sun: transition region 



1. INTRODUCTION 

The Sun displays a bewildering array of jet-like phe- 
nomena that can be observed at many different wave- 
lengths and in many different regions. At the limb, 
we find spicules, protrusions of cool gas that can reach 
heig hts of 6-10 Mm and have lifetimes of several minutes 
(e.g. lBeckersllT968l ). and that are observed in chromo- 
spheric lines such as Ha and Ca Ii H. There are also larger 
and longer-lasting jets such as surges and macrospicules. 
Recently, faster and shorter -lived jets, called type I I 
spicules, have been found ()De Pontieu et al.l l2007bf ). 
many lasting less than 100 s before fading from view. 
Hinode observations were instrumental in finding these, 
but have also shown a large am ount of other jet activ- 
ity in Ca H and ot her Hues (e.g. iKatsukawa et al.ll2007t 
IShibata et"al]|2007f ). 

On the disk, mottles show much of the same behaviour 
as spicules, and we can also observe shorter-lived dy- 
namic fibrils in the same lines. In higher-temperature 
ultraviolet lines such as C iv, Si iv, and O vi, we observe 
short-lived explosive events and longer-lasting blinkers. 
The former, in particular, have ve ry wide line prof iles 
indicating strong bidirectional jets (llnnes et al.lll997l ). 

There have been investigations into whether some 
of these terms simply describe different aspects of the 
same basic phenomena, e.g. spic u les and mottles 
(iGrossmann-Doerth fc Schmidtl 119921; pTsiropoula et al.l 
Il99 4f). blinkers and explosive j3vents ("Ch ae et al.ll200dl ). 
spicules and explosive events (iWilh elm 2000l). or spicules, 
blink ers, and explosive events (jMadiarska fc Dovld 
l2003l ). Different authors often conclude differently, and 
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the matter can not be considered settled. Part of the 
problem is that we usually lack cotemporal, cospatial 
observations across several wavelength bands, so that we 
most often cannot trace events across different tempera- 
tures and with enough spatial and temporal resolution to 
elucidate the relationship between the various features. 
Another complication is that we cannot observe the ex- 
act same phenomenon both on the disk and at the limb, 
for simple geometric reasons. Simulations can be very 
helpful here in letting us change our vantage point at 
will. 

The driving mechanism of the various jets is in gen- 
eral not well established. The main candidates have 
been waves (magnetoacoustic or Alfven) and magnetic 
reconnection. There are strong indication s that acoustic 
Shockwaves cause dynamic fibrils (Hansteen et al.l 120061 : 
iDe Pontieu et al.ll2007al: iHeggland et al.ll2007D and that 
reconnection cau ses explosive events (jinnes et all 119971 : 
IChae et all Il998[ ). but the case is more o pen for the 
other type s of jets. It has been noted by e.g . IChae etldl 
([1998), N ine et al.l (|2004[ ). and Dov le et alj (|2006D that 
many explosive events occur in bursts, often with inter- 
vals of 3-5 minutes. These periods correspond to the 
dominant wave modes produced by the solar granulation, 
suggesting that such waves may induce or modulate re- 
connection events. 

A rey i ew of various spicule models can be found in 
ISterlind (|2000f) . Most models so far are in ID, assum- 
ing a rigid magnetic flux tube either with or without 
expansion. Some use a piston driver, either at a given 
frequency or randomised, while others use a sudden in- 
crease in pressure and temperature as a trigger mecha- 
nism for setting up shocks. This pressure increase is of- 
ten as sum ed to be in the photos phere, but IShibata et al.l 
(|1982[ ) and lSterhng et all ()1993D have investigated the ef- 
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fects of energy input in the upper chromosphere as well, 
producing different types of jets that have similar prop- 
erties to spicules or surges. The source of the energy 
input is usually not specified, though reco nnection could 
be a n atural candidate, as speculated by [Sterling et al.l 

2D simulations containing an actual magnetic field (as 
opposed to treating it m erely as a given rigid flux tu be) 
have been performed bv iTakeuchi fc Shibatal (|2001al lH). 
They studied photospheric reconnection which they 
claim produces a large enough wave energy flux to drive 
spicules, but their actual simulation box does not extend 
past the lower chromosphere. 

Explosive e vents have been model le d in ID by 
lErdelvi et all ([199^ and iSarro et all (fTool . 2D 
simulations of explosive eve n ts have bee n per - 
formed by iKarpen et a l.l (I1995D llnnes fc Totf 



, jd] (fl995l) llnnes fc Tothl 

iRoussev et al.l pOOlallbU d). and IChen k Priesti 



1999 ). 



20061) 



Again, the ID simulations generally assume a rigid flux 
tube geometry, and treat the magnetic reconnection 
not explicitly, but as a sudden deposition of energy at 
a specified height. The 2D simulations often assume a 
simplified ma gnetic geometry with vertical antiparallel 
field lines. Ka rpen et al.l (|1995D use a more complicated 
magnetic field geometry with some similarities to the 
one used in this paper, but do not include radiative 
losses or heat conduction in their simulations. 

In this paper, we present the results of 2D simulations 
that include the effects of radiation and heat conduction, 
and that involve a complex magnetic field geometry in 
which reconnection events are induced by waves prop- 
agating upwards from the photosphere/convection zone. 
We do not attempt to match the results to specific obser- 
vations, but show the signatures of these events in several 
chromospheric, transition region and lower coronal spec- 
tral lines, and viewed both on the disk and at the limb 
of the Sun. 

2. SIMULATIONS 

The model a tmosphere is similar to the one used by 
iHeggland et~al (2007), but now the lower boundary of 
the domain reaches down to the upper photosphere, as 
opposed to the chromosphere. The simulation box ex- 
tends about 11 Mm in height, going from the photo- 
sphere through the chromosphere, transition region and 
lower corona. The lower boundary, z — 0, corresponds to 
a hei ght of 150 km in the VAL3C model (|Vernazza et al.l 
Il981| ). The simulation box contains 201 x 191 grid cells, 
using a uniform spacing of 50 km in the horizontal direc- 
tion and a non-uniform spacing in the vertical direction, 
starting at 16 km in the chromosphere and transition 
region and increasing exponentially in the corona, reach- 
ing 220 km at the upper boundary. In Figure [1] we have 
plotted the initial temperature distribution with mag- 
netic field lines and curves of equal plasma /3, the latter 
defined as the ratio between the thermal and the mag- 
netic pressure. 

The magnetic field configuration is a potential field, 
with a central flux tube which progressively widens with 
height and ends up dominating the field in the corona. 
However, there are also two inclined flux tubes, one on 
each side of the centre. These tubes both end up in 
magnetic null points, located at transition region heights 
{z = 1.75 Mm, X = 2.6 Mm and x = 7.4 Mm). The nuU 




Fig. 1. — Initial temperature structure, with overplotted /3 con- 
tours [white) and magnetic field lines (green). The red line between 
X = 4 and x = 5 Mm marks the location of the piston. The mag- 
netic null points show up as local maxima of the /3 at 2 = 1.75 Mm 
and X = {2.6,7.4} Mm. 
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Fig. 2. — Plot of the vertical velocity, set to saturate at 1 km s~ . 
Selected magnetic field lines in green, Alfven speed contours in 
km s~^ in white. The first fast mode wave reaches the transition 
region after 12 seconds, refracts away from the high Alfven speed 
in the centre, and perturbs the field in the vicinity of the magnetic 
null points, triggering reconnection. 

points can be seen in Figure |T] as (3 maxima that are very 
small in spatial extent. It is near these null points that 
the reconnection events we study in this paper occur. 
Finally, there are two outer flux tubes that open into the 
corona. 

It should be noted that this null point configuration 
is somewhat peculiar to the 2D geometry, and would be 
difficult to reproduce exactly in a 3D model, or indeed 
on the real Sun. However, it is our belief that the ob- 
servational signatures of reconnection events on the real 
Sun should not be strongly dependent on the specific 
geometry of the magnetic field, and that our simulated 
observations should be broadly similar to what we would 
expect to see when observing actual reconnection events. 
This is supported by the similarities we find with obser- 
vations of several types of events thought to be caused 
by reconnection. 

The co de we use is the sam e 3D MHD code that was 
used in.Heggland e t ahl (I2007T) . and is described in more 
detail by iHansteeni ( 20051 ). It should be noted that it 



contains a realistic radiative loss function based on coUi- 
sional excitation of hydrogen, carbon, oxygen, neon, and 
iron, as well as thermal c onduc tion along the magnetic 
field (see 'Hegg land et al.l (|2007l ) for the details). 

In the models presented here we operate with a mag- 
netic diffusivity rj many orders of magnitude larger than 
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that on the Sun, and dissipation starts at much smaller 
magnetic field gradients. The dissipated energy is 

Qjoulo = E ■ J, (1) 

where J = V x B is the current density and the resistive 
part of the electric field is given by 

E^. = [l{vi'^+Vi'') + I{vi''+Vi'')]j., (2) 
and similar for Ey and ■ The diffusivities are given by 
V^'^^P^iviCf+V2\u,\) (3) 

= p^t;3|V^.u|_, (4) 

where PrM is the magnetic Prandtl number, c/ is the 
fast mode speed, wi, V2 and W3 are dimensionless numbers 
of order unity and the other symbols retain their usual 
meanings. 

The upper coronal boundary is maintained at 1 MK, 
while radiation and heat conduction set the temperature 
structure in the rest of the domain. 

In the simulations, we investigate the effects of different 
locations and methods of driving, and how the produced 
waves can trigger reconnection. In one case, we use a 
monochromatic localised piston with 300 s period and 
1.1 km s~^ amplitude, driving the lower boundary in the 
vertical direction at 2 = Mm between a; = 4 Mm and 
X — h Mm (marked with a red line in Fig.[T]). In another, 
we use a driver taken from observations with Hinode/SP 
(Lites, private communication), producing waves at all 
locations along the lower boundary and in a more real- 
istic range of periods. The data were obtained on 2007 
October 22 and contain the line-of-sight velocities in a 
slice across a network element. 

We also experiment with putting the null points at 
different heights in the atmosphere. In one case, they 
are located in the transition region at a height of about 
1.75 Mm above the lower boundary, as in Figure[Tl in an- 
other, they are in the upper chromosphere, some 450 km 
lower. 

3. ANALYSIS 

3.1. Transition region null points, piston driver 

3.1.1. General description of events 

We first look at the case with null points in the tran- 
sition region and the localised piston driver as described 
above. Due to the varying magnetic inclination in the re- 
gion of the piston, it will generate a mixture of fast and 
slow mode waves. In low-/? regions, which include the 
piston location as well as most of the simulation box, the 
fast modes can propagate everywhere, though they tend 
to be refracted into regions of low Alfven velocity (e.g. 
lOsterb rock 1961), such as null points. The slow modes 
are mainly restricted to propagation along the magnetic 
field. 

The waves disturb the field as they travel, and once the 
disturbance reaches the null points, reconnection is trig- 
gered as lines of opposite polarity are pushed together 
and meet. This process releases a significant amount 
of energy in a short time and leads to very rapid heat- 
ing of the plasma; in the stronger events, the plasma 



around the null points reaches coronal temperatures of 
roughly 1 MK. In addition, a bidirectional jet is formed, 
as plasma is rapidly accelerated away from the recon- 
necting region. 

Since the lines of opposing polarity are located quite 
close to each other around the null point, only a small 
disturbance is required to start off the reconnection pro- 
cess the first time. This happens soon after the start of 
the simulation, as fast mode waves generated by the pis- 
ton in the central flux tube propagate upwards and reach 
the transition region in as little as 10-12 seconds. There, 
they are refracted away from the high Alfven speed re- 
gion in the centre of the box (see Fig. [2|) and move in- 
stead towards the sides, where they reach the magnetic 
null points and push the field lines together, triggering 
reconnection. 

The piston, driving generally along the field in a mag- 
netically dominated (i.e. low-/3) plasma, is more efficient 
at generating slow mode than fast mode waves, so these 
initial fast disturbances have low amplitude — around 
1 km s^^ just before reaching the transition region — 
but this is enough to push the field lines above and be- 
low the null point into each other and trigger a weak re- 
connection event. This leads to moderate heating of the 
reconnecting region and the formation of a bidirectional 
jet, with high velocities upwards to the left and down- 
wards to the right. The hot jet going left then pushes 
cooler material sideways and upwards, propelling a small, 
spicule-like protrusion of cooler chromospheric material 
into the corona (Fig. [Sj left panel) . 

After some 160 seconds, or half a driver period later, 
the piston driven disturbance changes sign, and starts 
pulling the field lines above and below the null point away 
from each other. Instead, the lines to the left and right 
of the null point are pushed together, and reconnection 
occurs again, this time in a direction perpendicular to the 
original reconnection event. Another bidirectional jet is 
formed, this time being mostly vertical, and the region 
around the null point (now embedded entirely in cooler 
chromospheric gas) is heated to coronal temperatures. 
This heated region spreads out along the jets and the 
magnetic field, and creates a butterfly-shaped region of 
heated plasma (Figure [3l left centre panel) . 

The orientation of the upwards-moving jet changes as 
the reconnection continues, and it rotates from being di- 
rected upwards to the right (mainly across the field lines) 
to being directed somewhat to the left (mainly along the 
field lines) . The latter orientation allows the gas to move 
much more freely, and the hot plasma is soon accelerated 
upwards along the field lines with great speed (Fig. [3l 
rigtit center panel). As it passes the now mainly side- 
ways moving cool jet generated by the first reconnec- 
tion, it pushes some of that cool gas upwards as well and 
greatly extends the length of the cool jet, which reaches 
a maximum height of z = 8 Mm. 

This pattern then repeats itself, with reconnection be- 
ing triggered every 150 s (i.e. every half period), alter- 
nating between vertical and horizontal orientation, and 
propelling both cool and hot jets upwards. The later 
events are not as strong as the first one, and are quite 
similar to each other. An example of a later horizontally 
oriented reconnection jet is shown in the right panel of 
Figure [3l 
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Fig. 3. — Plots of the plasma temperature (greyscale) with magnetic field lines (green) and velocity vectors (red) at four different times 
in the model with transition region null points and piston driver. The plots zoom in on the area around the left-hand null point. At 74 s 
{left panel), a horizontal bidirectional jet is pushing cool material upwards to the left. At 206 s, reconnection has heated a significant area 
and propels a jet mainly across the field lines. At 236 s, the jet aligns more with the field and hot material is ejected upwards. The final 
panel shows a later horizontally oriented reconnection after 360 s. The field lines in all panels are drawn to show the geometry rather than 
the strength of the magnetic field. 
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Fig. 4. — Plots of the Ca signal (top left), density (top right), tem- 
perature with magnetic field lines (bottom left), and field-aligned 
velocity (bottom right) in the case with transition region null points 
and piston driver. The elapsed time is 226 s. Reconnection is pro- 
ducing a butterfiy-shaped, vertically oriented heated region and jet 
on the left, and a hot TR jet on the right. This figure is available 
as a movie in the electronic version of the Journal. 

3.1.2. Synthesised limb observations 

In order to get an impression of how these jets would 
look if observed at the solar limb, we have used the n on- 
LTE radiative transfer code MULTI (|Carlssonlll9"86h to 
calculate the population densities of the different excita- 
tion levels of Ca ii ions in our simulations. In the side 
views, we have chosen to concentrate on the Ca ii H line 
at 3968 A. In a 2D model, it is not possible to do a full 
radiative transfer calculation along a line perpendicular 
to the plane of the computational box, but the popula- 
tion density of the upper level of this transition serves as 
a rough approximation of the actual intensity. 

In Figure IH we show plots of the logarithm of the cal- 
culated Ca II H upper level density {upper left, henceforth 
and in the figure captions referred to as the Ca signal), 
the logarithm of the mass density from the MHD simu- 
lation data (upper right), the plasma temperature with 
overplottcd magnetic field lines [lower left), and the ve- 
locity along the magnetic field [lower right) after 226 s. 
At this time, reconnection with vertically oriented jets 



Fig. 5. — Plots of the Ca signal (top left), density (top right), tem- 
perature with magnetic field lines (bottom left), and field-aligned 
velocity (bottom right) in the case with transition region null points 
and piston driver. The elapsed time is 364 s. The spicule-like cool 
jet on the left has reached its greatest length, extending nearly 
8 Mm above the photosphere. 

is happening near the left-hand null point, while a hori- 
zontally oriented event is happening near the right-hand 
null point — though the jet on the right-hand side of 
it is field-aligned and has turned upwards quite sharply. 
The cooler jets show up in Ca ii as rather thin, elon- 
gated features, similar to spicules or fibrils. In contrast, 
the hotter jets, and the reconnection regions themselves, 
show no Ca ii signal. The hot jet on the right does have 
a density enhancement [upper right panel), but it and 
the regions around the null points are heated to temper- 
atures of about 1 MK, far above the temperature where 
calcium gets multiply ionised, which explains the lack of 
Ca II emission. The high velocities of the reconnection 
jets are clearly visible in the lower right panel. 

The cooler jets are seen to rise and fall; in Figure El 
we show the situation after 364 s, when the cool jet seen 
on the left in Figure |3] has reached its maximum extent, 
after being pushed upwards by the hot jet from the re- 
connection event at 200-240 s (ref. Fig. [3]). It reaches a 
height of nearly 8 Mm. 

In FigurelH we see a later, thinner cool jet on the left. 
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Fig. 6. — Plots of the Ca signal (top ie/t), density (top right), tem- 
perature with magnetic field lines {bottom left), and field-aligned 
velocity (bottom right) in the case with transition region null points 
and piston driver. The elapsed time is 772 s. A thin, low-density 
cool jet can be seen on the left in the temperature plot. 

which shows up only weakly in Ca ii. The density plot 
(upper right) tells us why: its density appears to be too 
low, even though it is in about the right temperature 
range. 

The typical lifetimes of the cool jets are 200-300 s, and 
they reach maximum lengths of 6-8 Mm above the photo- 
sphere. Their thickness varies a lot with time, but is typ- 
ically slightly less than 1 Mm. These fig ures are roughl y 
within the ranges reported for spicules (jBeckerslll968f ). 
However, it is important to look at the whole range of 
spectral, spatial and temporal data available before try- 
ing to establish correspondence between simulations and 
observations. 

3.1.3. Synthesised disk observations 

In producing synthesised disk observations, we look 
into the simulation box directly from above, in the plane 
of the box. A proper radiative transfer treatment is then 
possible, and has been carried out for the Ca ii IR line 
at 8542 A, widely used in spectroscopic studies. The 
calculations have been made column by column, i.e. ne- 
glecting any radiative interaction in the x-direction. A 
simplified treatment based only on coUisional excitation 
has been carried out for the optically thin lines of C iv 
(1548 A), O VI (1032 A), and Fe xii (195 A). These lat- 
ter lines will typically be formed in the transition region 
or lower corona, whereas the Ca IR line is formed in the 
chromosphere. In the treatment of the Fe line, which 
has a wavelength shorter than the Lyman threshold at 
912 A, we have included the effects of absorption by neu- 
tral hydrogen and neutral and singly ionised helium. 

A-t plots of these lines, looking directly down on the 
centre of the reconnection region a,t x = 2.6 Mm, are 
shown in Figure [T] The C, O, and Fe lines have had 
their resolution and cadence downgraded to typical in- 
strumental values of 650 km (about 1") and 20 s, and 
are plotted using logarithmic scaling to bring out details 
in the weaker later events. Ca is plotted at 150 km res- 
olution and 2 s cadence, using linear scaling. 
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Fig. 7. — Spectra of the high null point piston driver case, cen- 
tered at X = 2.6 Mm and with blueshifts shown as negative velocity. 
The C, O and Fe lines use logarithmic scaling in order to bring out 
the details of the much weaker later events. The difference between 
white (maximum intensity) and black (intensity cutoff) is a factor 
10'*. The UV lines experience strong brightening and extensive 
Doppler shifts when the reconnection events happen. The profiles 
are asymmetric, having higher blueshifts than redshifts, but the 
intensity maximum is slightly redshifted. The Ca spectrum uses 
linear scaling and a much smaller Doppler shift scale. 

In the transition region lines, there is a sudden onset 
of a very bright, very broad feature, with upflows of 80- 
90 km s~^ and downflows of 60-80 km s~^. This happens 
first in C iv, the lowest temperature line, at 160 s, then 
in O VI at 180 s (with a weaker signal at 160 s). This fea- 
ture is caused by the sudden heating and the powerful 
bidirectional jet generated by the reconnection. These 
spectra sho w strong similarities to observations o f explo- 
sive events (jlnnes et all 119971: IChae et all |2000D which 
are believed to be caused by reconnection. iDerg ()1994l ) 
finds that explosive events have typical lifetimes of 60 s, 
spatial extents of about 1.5 Mm, and typical Doppler 
shifts of 100 km s^^, though both higher and lower dura- 
tions and Doppler shifts have been reported. The figures 
match up well with our results; the main parts of the 
brightcnings fade away after 60 s and the spatial extent 
is roughly in the quoted range (see Figs. [3] and |4]). 

The profiles are asymmetric, with larger blueshifts 
than redshifts, and in particular the C and O lines have a 
short-lived blueshift of 140-150 km s'^ at 220 s. On the 
other hand, the intensity maximum is redshifted. This 
behaviour can be readily understood in terms of the den- 
sity stratification. The reconnection outfiow jets are ex- 
pected to propagate outwards at the Alfven speed, which 
is inversely proportional to the square root of the den- 
sity. Hence, the jets propagating upwards in a stratified 
medium should be faster than those propagating down- 
wards, leading to an asymmetry like the one we find. The 
blueshift maximum between 220 and 240 s happens when 
the heated jet, which is initially embedded in material at 
chromospheric temperatures (Fig. [4|) , enters the corona 
where the Alfven speed increases markedly. Meanwhile, 
the compression effect of the jet should be stronger for 
the part propagating downwards into denser material. 
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leading to higher emission at shght redshifts. 

The temperature in the reconnection region gets high 
enough to produce significant emission also in the Fe xii 
line, but it is quite strongly affected by neutral hydro- 
gen absorption. The first powerful event is completely 
absorbed because the reconnection region at that time 
is located under a small jet with at least 1 Mm of chro- 
mospheric material above it. The first horizontal jet re- 
connection event shows up weakly at around 400 s, but 
strong blue- and redshifts are not seen until the vertical 
jet at 500 s. The intensity maximum is in a horizontal 
jet at around 680 s, at which time there is hardly any 
chromospheric material above the reconnection region. 

The chromospheric Ca ii IR line shows signs of activity 
before any of the higher temperature lines. A relatively 
strong upflow (15 km s^^) is caused after less than 100 s 
by the first weak reconnection episode. This upfiow is 
soon reversed, and after 150 s the first major reconnec- 
tion happens, changing the line centre from a redshift 
of 20 km s~^ to almost instantly, while also yielding 
a significant blueshifted component with speeds up to 
30 km s-i. 

After this violent first event, several periodic weaker 
reconnection events follow. In these, the line only un- 
dergoes significant redshifts, with very little evidence 
of blueshifts. This indicates that the line's main for- 
mation height is below the reconnection region. The 
other lines show periodic weaker brightenings, though 
still with quite high velocity. In every other event, the 
main jets are oriented horizontally, leading to smaller 
Dopplcr shifts. The intensity decrease is mainly due to 
reduced density in the reconnection region; the reconnec- 
tion jets cause a large outflow of material, and because 
the inflow (amplified by the waves) soon triggers another 
reconnection event, the area never recovers its initial den- 
sity. 

3.2. Chromospheric null points, piston driver 
3.2.1. General description 

These simulations use a piston driver of the same pe- 
riod (300 s) and amplitude (1.1 km s^"'^) as in the previ- 
ous section. However, here the magnetic null points are 
located around 450 km lower in the atmosphere, putting 
them in the upper chromosphere rather than in the lower 
transition region. The basic geometry of the magnetic 
field is not changed, as we can see in Figure [8l the only 
difference is a translation of the whole geometry to lower 
heights. 

3.2.2. Synthesised limb observations 

Figure [8] shows the simulation, viewed from the side, at 
the time when the first and most powerful reconnection 
event has reached its maximum, at around 260 s elapsed 
time. Once more, we have plotted the Ca H signal {up- 
per left), the logarithm of the mass density (upper right), 
the temperature [lower left) and the field-aligned velocity 
{lower right). We clearly see the butterfly-shaped region 
of higher temperature {lower left) and the fast bidirec- 
tional jet {lower right), although the latter is now going 
upwards much faster than downwards due to the plasma 
being able to more easily propagate along the magnetic 
field. The beginnings of a jet can be seen in Ca H {upper 
left) , although the heated region is too hot (and to some 
extent too evacuated) to show any Ca signal. 
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Fig. 8. — Plots of the Ca signal {top left), density {top right), tem- 
perature with magnetic field lines {bottom left), and field-aligned 
velocity {bottom right) in the case with chromospheric null points 
and piston driver. The elapsed time is 260 s. The outflow speed 
from the reconnection region on the left is at its maximum. This 
flgure is available as a movie in the electronic version of the Journal. 
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Fig. 9. — Plots of the Ca signal {top left), density {top right), tem- 
perature with magnetic field lines {bottom left), and field-aligned 
velocity {bottom right) in the case with chromospheric null points 
and piston driver. The elapsed time is 354 s. The extent of the 
cool jet is at its maximum, 4 Mm above the photosphere. 

In Figure [H we show the situation as the cool jet 
reaches its maximum extent, 90 seconds later. As could 
be expected, the velocity signal is weak at this point, and 
the region around the null point is no longer heated to 
coronal temperatures. The jet shows up clearly in Ca, 
but reaches a maximum height of only about 4.5 Mm, 
as compared to 8 Mm when the null points are located 
in the transition region. The general pattern described 
above repeats itself periodically every 300 s, with pro- 
gressively slower and shorter jets. 

Overall, this case is quite similar to the case with 
higher null points described in the previous section. The 
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Fig. 10. — Spectra of the chromospheric null point piston driver 
case, centered at a; = 2.6 Mm and using logarithmic scaling for 
the UV lines and linear scaling for Ca. C and O show largely 
symmetric profiles with maxima on each reconnection event with 
vertical jets. Fe shows only the background emission. Ca has a 
complex line profile with multiple components. 
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Fig. 11. — Ca spectrum with doppler shifts (grey scale) in the 
chromospheric null point piston driver case, with velocities at z = 
0.93 Mm {solid line) and z = 1.37 Mm (dashed line) superimposed. 
The former fits well with the line centre position, while the latter 
is a reasonable match for the blueshifted spikes and their following 
redshifts. 

jets are slower and shorter because of the higher density 
surrounding the nuU points, resulting in lower accelera- 
tion for similar force. 

3.2.3. Synthesised disk observations 

Figure [To] shows A-t diagrams of C iv, O vi, Fe xii and 
Ca II IR as they would appear looking straight down at 
the centre of the reconnection region {x — 2.6 Mm), as 
for the high null point case in section 13.11 The resolu- 
tion is also the same as used previously. In C iv and 
O VI, we get bidirectional signals spaced 300 s apart. 



with Doppler shifts of 40-60 km s^^, but little other sig- 
nal. In Fe XII, the events give no visible signal, and we 
just see the background emission. The reason for the low 
Fe emission is a combination of the relatively lower tem- 
perature reached in the region heated by reconnection 
and greater absorption by neutral hydrogen due to the 
deeper location of the reconnection region. The UV lines 
are in general more symmetric in this simulation, because 
the hot jet does not reach the perturbed transition region 
and corona before being cooled. 

The Ca line, though, appears highly complex in this 
simulation. It exhibits both absorption and emission, as 
well as strong Doppler shifts into both red and blue. In 
observations, such a spectrum would be difficult to make 
sense of. However, we have the full simulation data to 
work with, including the velocity at all heights. In fact, 
at least two different velocity components can be identi- 
fied as being responsible for the appearance of the spec- 
trum, as shown in Figure 1111 The line centre closely 
follows the vertical velocity at a height of 0.92 Mm, be- 
low the magnetic null point, as marked with a solid white 
line in the figure. The strongest blueshifted signal, which 
later passes through the line centre and becomes red- 
shifted, corresponds reasonably well to the vertical ve- 
locity at a height of 1.37 Mm [dashed line), just above 
the null point, which is located at 1.30 Mm. These two 
velocity components are in counterphase as a result of 
the bidirectional jet produced by the reconnection. 

This nicely illustrates the fact that spectral lines are 
formed over a range of heights, rather than at one specific 
height, and that care must be taken when using them as 
probes of atmospheric conditions. But it also illustrates 
that such a multi-component signal could be used as evi- 
dence of reconnection via the characteristic bidirectional 
velocity pattern, as long as the reconnection happens in 
the height range where the line is formed. By contrast, 
in our simulation with null points in the transition re- 
gion, we observe only the downflows because the line is 
primarily formed below the null points and reconnection 
region. 

3.3. Transition region null points, Hinode driver 
3.3.1. General description 

In this simulation, we go back to the magnetic configu- 
ration of the first simulation studied (section [3.ip . where 
the null points are located in the transition region. How- 
ever, instead of using a localised, monochromatic piston 
driver, we use velocity data obtained with Hinode/SP (as 
described in the introduction) to drive the whole lower 
boundary. As before, we will not attempt to match the 
simulation to specific observations; the data simply serve 
as a more random and realistic velocity driver with a 
wider range of frequencies than our monochromatic pis- 
ton. 

3.3.2. Synthesised limb observations 

Using the Hinode driver changes a number of things. 
In general, the velocity field in the atmosphere becomes 
much more complex. On the left-hand side, the first 
reconnection has its jet oriented vertically, as opposed 
to the horizontal orientation we get in the runs with a 
piston driver. This is a result of the different phase of the 
first waves to reach the null point. In this case, as the 
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Fig. 12. — Plots of the Ca signal {top left), density (top right), 
temperature with magnetic field lines {bottom left), and field- 
aligned velocity {bottom right) in the case with transition region 
null points and Hinode driver. The elapsed time is 184 s. A strong 
and very fast hot jet is being produced by the reconnection at the 
right-hand null point. This figure is available as a movie in the 
electronic version of the Journal. 
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Fig. 13. — /3 plot showing the height change of the right-hand 
null point with time when it starts out in the transition region and 
the Hinode driver is applied. 

reconnection happens in the transition region, only the 
hotter material above the null point is ejected upwards, 
and there is little cool material above the null point to 
be accelerated by later reconnection events. As a result, 
only few and short cool jets are formed at all. 

At the right-hand side null point, no significant recon- 
nection happens until 130 s, but then a quite major event 
occurs, with its jet oriented vertically. This leads to the 
ejection of a powerful hot jet, as seen in Figure [T2l The 
jet is too hot to show up in Ca. 

This latter reconnection event is actually quite long- 
lasting, the flow direction only being reversed at 290 s, 
and then only for 90 s before another 170 s period of verti- 
cally oriented reconnection. The sustained event releases 
quite a bit of energy from the magnetic field and changes 
its structure fairly significantly; among other things, the 
null point is moved down into the chromosphere by 3- 
400 km (Fig. [13)) . Later jets are fairly weak, the cool 



Fig. 14. — Plots of the Ca signal {top left), density {top right), 
temperature with magnetic field lines {bottom left), and field- 
aligned velocity {bottom right) in the case with transition region 
null points and Hinode driver. The elapsed time is 326 s. This 
figure shows the longest of the cool jets that appear in this simu- 
lation. 
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Fig. 15. — Spectra of the high null point Hinode driver case, 
centered at x = 7.6 Mm. The UV lines use logarithmic scaling and 
the Ca line uses linear scaling. The C and O lines show mainly 
redshifts because the matter moving upwards is heated out of their 
passbands. Fe shows not only the signal of the heated area around 
the reconnection point, but also of the fast hot jet to the side of it, 
shown in more detail in Fig. 1161 

ones reaching only around 1 Mm above the transition 
region (Fig. [Tl)) . 

The long duration of the event is likely a result of the 
specifics of the driver. On the right-hand side of the box, 
the input velocity only rarely changes sign, leading to a 
fairly continuous flow in one direction, rather than the 
periodic reversals of the piston driver. 

3.3.3. Synthesised disk observations 

Because the null point in this simulation displays sig- 
nificant movement, both horizontally and vertically, and 
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Fig. 16. — Fe XII spectrum of the high null point Hinode driver 
case, centered at x = 8.45 Mm and using logarithmic scaling. This 
spectrum shows the very fast hot jet and its considerable blueshift 
of 160 km s-i. 
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Fig. 17. — Ca II IR spectrum of the high null point Hinode driver 
case, centered a t x = 7.9 Mm. Weak blueshiftcd excursions similar 
to those in Fig. llOl can be seen as a result of the null point having 
moved down after extended reconnection. 

because the jets tend to follow the magnetic field, which 
is slightly inclined even in the lower corona, we can find 
interesting phenomena at several different locations when 
looking down as we would do if the events happened on 
the solar disk. We calculate the same C, O, Fe, and Ca 
lines as above, and use the same spatial resolution and 
exposure times. 

In Figure [HI we show A-t plots of the different spec- 
tra centered at x = 7.6 Mm. This is the place where 
the transition region and coronal lines have their high- 
est intensities, looking straight down on the reconnec- 
tion region. The C iv line shows little bidirectionality 
in this case, having its maximum at redshifts of around 
25 km s~^. The O vi line is bidirectional, but has much 



higher intensity in the red. The Fe xii line is more ex- 
tended, but also has its maximum at the same redshift. 
This line not only includes emission from the region im- 
mediately surrounding the null point, but also from the 
upwards propagating fast hot jet between a; = 8 and 
a; = 9 Mm. A spectrum centered &t x = 8.45 Mm show- 
ing this jet in more detail is shown in Figure [TBI We can 
see the strong signature of the fast hot jet and its decel- 
eration between 160 and 340 s. It reaches an impressive 
blueshift of roughly 160 km s~^, matching the maximum 
field-aligned velocity from the simulation data. 

The Ca il IR line shows only redshifts for this event 
(130-200 s), some of them rather powerful (up to 
20 km s~^), and exhibits strong central reversal at 
around 140 s, when the plasma is significantly heated by 
the reconnection. The later Ca spectrum is more com- 
plex, likely showing several overlapping velocity compo- 
nents as in the low null point piston driver case in sec- 
tion 13.21 as well as the effects of the movement of the 
null point itself. Of course, the velocity input from the 
Hinode driver includes a wide spectrum of frequencies 
rather than the single one of the piston, which serves to 
further complicate the simulated Ca spectrum. 

The effects of having the null point move several hun- 
dred km downwards are visible in the Ca signal in some 
locations. In Figure 1171 we show a portion of the spec- 
trum at X = 7.9 Mm. Between 680 and 780 s, there is 
a clear blueshifted excursion, similar to what we saw in 
the spectra of the simulations with chromospheric null 
points (Figs . [TOl and [TT|) . By the time of this event, the 
right-hand null point in this simulation is at a height of 
1.4 Mm (Fig. [T3[) . very close to the height of the null 
points in the other simulation. 

3.4. Chromospheric null points, Hinode driver 
3.4.1. General description 

In this simulation, we again use the magnetic configu- 
ration that puts the magnetic null points in the chromo- 
sphere, but use the same Hinode driver as in the previous 
section. Once again this leads to extensive reconnection 
near the right-hand null point, but because that is now 
embedded in chromospheric plasma, it produces a major 
cool jet rather than the hot one produced when the null 
points are in the transition region. Again the reconnec- 
tion changes the field and pushes the null point lower in 
the atmosphere after about 500 s fFig. [T5)) . Later jets are 
fairly weak and do not extend far above the transition re- 
gion; neither do the jets produced by the weaker events 
at the left null point. Some jets are also produced by 
shock waves propagating from the lower boundary, espe- 
cially in the magnetic channel connecting to a; = 2 Mm 
— the shock fronts are clearly visible in movies of the 
field-aligned velocity. This pr ocess is similar to th e one 
that produces dynamic fibrils ([Hansteen et al.ir2006f ). but 
the jets produced this way tend to reach heights of no 
more than 0.5 Mm above the transition region in this 
simulation. 

3.4.2. Synthesised limb observations 

In Figure [TH we show the Ca H signal, density, tem- 
perature and field-aligned velocity at 216 s elapsed time, 
when the region heated by the first powerful reconnection 
event has reached its maximum extent. We see that the 
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Fig. 18. — /3 plot showing the height change of the right-hand 
null point with time when it starts out in the chromosphere and 
the Hinode driver is applied. 
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Fig. 20. — Plots of the Ca signal {top left), density (top right), 
temperature with magnetic field lines {bottom left), and field- 
aligned velocity {bottom right) in the case with chromospheric null 
points and Hinode driver. The elapsed time is 350 s. The cool jet 
has reached its maximum extent, 8 Mm above the photosphere. 

and a duration P of roughly 430 s. Altho ugh these values 
fall outside the range studied by Heggla nd et al.l (|2007D 
in their study of shock wave-driven dynamic fibrils — in 
particular, the maximum velocity is much higher than 
would reasonably develop through steepening of acoustic 
waves alone — they are a pretty good fit to the formula 
predicted and found in that paper, namely 



d = 



P/2 ' 



(5) 



K [Mm] 



X [Mm] 



not too surprisingly since the jet is driven by a shock, 
albeit one generated by reconnection rather than directly 
by photospheric motions. 



Fig. 19. — Plots of the Ca signal {top left), density (top right), 
temperature with magnetic field lines (bottom left), and field- 
aligned velocity {bottom right) in the case with chromospheric null 
points and Hinode driver. The elapsed time is 216 s. A powerful 
reconnection event is happening at the right-hand null point, heat- 
ing a large area and propelling a cool jet upwards at great speed. 
This figure is available as a movie in the electronic version of the 
Journal. 

upper part of that region follows the field lines upwards; 
the total length of the heated region is more than 2 Mm. 
We also see the high velocity of the jet. Note that the 
heated region gets too hot to show up in the Ca signal. 

In Figure [201 we show the situation at 350 s, when the 
cool jet has reached its maximum extent. As we see, it 
reaches a height of nearly 8 Mm above the lower bound- 
ary, or 6 Mm above the transition region. This cool jet 
shows up quite strongly in Ca, and has a noticeable den- 
sity enhancement compared to the surrounding corona. 
The velocity, as one would expect, is close to zero at this 
time. 

This jet does not form a perfect parabola in a z-t dia- 
gram (it ascends slightly faster than it descends), but the 
best p arabolic fit, using the method of lDe Pontieu etldl 
(|2007a) and Heggland et all (|2007[ ). yields a maximum 
velocity Vmax of 52 km s""'^, a deceleration d of 250 m s~^. 



3.4.3. Synthesised disk observations 

In Figure!^ we show the simulated C iv, O vi, Fe Xil, 
and Ca ii IR spectra looking down on a: = 7.65 Mm. As 
in the other cases, the TR lines are plotted at 650 km 
and 20 s resolution, and the Ca spectrum at 150 km and 
2 s. 

C, O, and Fe all have huge maxima in connection with 
the first reconnection event, centered at 180-200 s — even 
with a lower intensity cutoff of lO"'* times the maximum 
intensity, as in all the UV spectral plots, only Fe shows 
any part of the background signal. The event shows up as 
a clear bidirectional jet in all three, with 50-70 km s~^ ve- 
locities; even up to 90 km s~^ in O VI. C and O have some 
of the same asymmetry as in the TR null point piston 
driver case f section l3.1|) . with slightly higher blueshifts 
than redshifts and slightly redshifted intensity maxima. 

In Ca, we get strong emission and a redshift of 
17 km for this event. Later on, the spectrum gets 
more complex, as usual with the Hinode driver. How- 
ever, it is notable that after the null point moves down at 
around 500 s, the Ca spectrum mainly shows blueshifts, 
since the primary formation height is now above the re- 
connection region. This makes a nice contrast with the 
mainly redshifted spectrum of the high null point pis- 
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Fig. 21. — Spectra of the chromospheric null point Hinode driver 
case, centered at x = 7.65 Mm. The UV lines use logarithmic 
scaling and the Ca line linear scaling. All lines have huge maxima 
when the first major reconnection event happens. Fe is largely 
symmetric while the C and O lines have asymmetries similar to 
those in Fig. [7] Ca, formed at the lowest height, is redshifted, but 
the later signal is shifted more towards the blue after the null point 
moves to lower heights. 
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Fig. 22. — Energy release rate in the form of joule heating in the 
four different cases: transition region null point with piston driver 
(black), chromospheric null point with piston driver (green), transi- 
tion region null point with Hinode driver (blue), and chromospheric 
null point with Hinode driver (red). 

ton driver case (Fig. [7]), where the nuU point is located 
almost 900 km higher up. 

3.5. Energy release estimates 

In Figure [2H we have plotted the energy release rate 
in the form of joule heating /a, where j is the current 
density and a the conductivity) in our four simulations. 
The energy release is calculated over the reconnection 
region, in this case defined as the region around the null 
point where the joule heating is above 5% of the value at 
the null point. The results are not very sensitive to the 
particular value of the threshold, and changing it from 
1% to 10% changes the calculated energy release by less 
than 10%. In the piston driver cases, we include the 



TABLE 1 

Energy release per reconnection event 



First event Strongest later event 

Case Energy [J] Duration [s] Energy [J] Duration [s] 

TR np/piston 1.23 X 10^* 280 5.40 X IQi'' 280 

Ch. np/piston 5.48 X lO^^^ 280 4.19 X IQi'' 280 

TR np/Hinode 5.01 X lO^^^ 240 7.15 X IQi'' 260 

Ch. np/Hinodc 3.31 X 10^^ 100 7.18 X 10^^ 280 

region around the left-hand null point, while we use the 
right-hand one in the Hinode driver cases. 

As we would expect from the side views and spectra, 
the greatest peaks in the energy release occur in the tran- 
sition region null point piston driver case {black line) 
and the chromospheric null point Hinode driver case {red 
line). The two piston driver cases are more regular, with 
double peaks (horizontal, then vertical jets) appearing 
every 300 s. The Hinode driver cases have second peaks 
appearing between 600 and 800 s, and the energy release 
is actually greater in these events than in the first one 
at around 200 s, but the later events show only weak ob- 
servational signatures. This is because the jets of these 
events are horizontally oriented, giving no Doppler shifts 
in the simulated spectra and not very tall jets, and also 
because the energy release happens over a longer period 
of time than at the first peak, at least in the case with 
chromospheric null points. 

The energy release rate is given in J m^^ s~^ because 
the total energy release will be dependent on the size of 
the reconnection region in the third dimension. In or- 
der to calculate the total energy release, we must assume 
an extent, and this value will necessarily be somewhat 
arbitrary. A value of 0.5 Mm has been assumed in Ta- 
ble [TJ The released energy per event is comparable to 
the generally assumed values for nanoflares, though not 
all the energy released in these events will be in the form 
of joule heating. 

4. DISCUSSION AND SUMMARY 

These simulations have been able to produce a number 
of jet phenomena through the same basic mechanism of 
wave-induced, or in the Hinode driver cases also convec- 
tive flow-induced, magnetic reconnection. In particular, 
they appear to be the first 2D simulations to show the 
formation of spiculc-like jets as a result of reconnection, 
while including most important physics; the main omis- 
sion is time-dependent ionisation. The simulated spicules 
match observed lifetimes and lengths. 

In this model of spicule formation, we also expect a 
brightening and broadening of UV lines at the footpoint 
of the spicule at the start of its lifetime; the heated 
plasma causing this is in general cooled down on a shorter 
timescale than the lifetime of the spicule. This UV bright 
point is a property our model shares with most ear- 
lier models that assume a sudden energy/pressure de- 
position as the source of spicules, rather than a veloc- 
ity perturbation at photospheric heights. Although UV 
bright points are frequently obser ved in connection with 
spicules, some observations (e.g. ISuematsu et al.l [l995[ ) 
indicate that this brightening happens after the ascend- 
ing phase of the spicule, rather than at the beginning. 
This timing problem is not resolved by our model. 

The UV brightenings that are produced in the recon- 
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nection regions have many of the same properties as ob- 
served explosive events, which could point towards a pos- 
sible connecti on between spic ules and explosive events, as 
suggested bv IWilhelml ()2000D . One of the plasma accel- 
eratio n me chanisms suggested in that paper, a slingshot 
effect (IWiIhelm.2000. Fig. 8) may be happening at times 
in our model; see for example the left centre panel of Fig- 
ure [3l where the mass flow is perpendicular to the field 
lines and chromosphcric matter is propelled upwards by 
several hundred km. On the other hand, flow along the 
field proves much more effective in lifting chromosphcric 
material in our simulations (right centre panel of Fig. [3] 
and lower two panels of Fig.|l]); different field geometries 
could give different results. 

Although our UV brightenings match many of the 
properties of explosive events, there are also some differ- 
ences. Our events are generally very bright, with max- 
imum intensities 50-100 times higher than the "quiet" 
profiles — in some cases, even 1000 times. Although 
some explosive events have strong brightenings, the ma- 
jority are not particularly bright Inncs (2001) quotes 
typical brightening factors of 2-5. Also, we get significant 
signal in the lower coronal Fe Xii 195 A lin e, while explo- 
sive ev ents rarely show up in coronal lines. lErdelvi et al.l 
()1999l ) show one example of an explosive event brighten- 
ing m TRACE 171 A, but since it does not show up in 
the lower temperature Mg X 625 A line, they interpret it 
as a brightening of a transition region O VI line within 
the passband rather than the coronal Fe X line. 

In one of our cases (Fig. [T0|) . we also find only very 
weak Fe xii emission. Since the heating is quite strongly 
localised to the reconnection region, it becomes a blob of 
coronal plasma surrounded by low-temperature chromo- 
sphcric material. In the low-temperature region, neutral 
hydrogen is present in significant proportions, and this 
leads to strong Lyman absorption of all lines with shorter 
wavelengths than 912 A. This, combined with the heat- 
ing events themselves being relatively weak, explains the 
lack of Fe emission in this case, and could also do the 
same for observed explosive events if in fact they occur 
relatively deep in the chromosphere. However, the C 
and O lines will be unaffected by this absorption, and 
are still stronger than observed. Also, the high veloci- 
ties observed in many explosive events indicate a quite 
high formation height. A more likely explanation for 
the discrepancy may therefore be that our reconnection 
events are simply more energetic than those responsible 
for most explosive events. The temperatures reached are 
high enough (> 1 MK) that they might even have an 
X-ray signal. 
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It should be noted that even the very strong first event 
in the TR null point piston driver case is completely ab- 
sorbed by the neutral hydrogen above it. Therefore, the 
absence of observed signal in high-temperature lines does 
not necessarily mean that the plasma is not heated to 
such temperatures, as the emission could be absorbed by 
overlying cool matter. 

Finally, wave-induced reconnection offers a natu- 
ral explanation for the repetitive be haviour sornetimes 
observed in exp losive events (e.g. IChae et all 119981 : 
iNing et al]|2004l ). If these events are caused by recon- 
nection induced by waves, we would naturally expect 
repetition at time scales of 3-5 minutes, corresponding 
to the dominant periods in the solar chromosphere and 
photosphere. 

Although we find a number of similarities between our 
simulation results and observations, the complexity of 
the signals in our synthetic spectra and images shows 
that it is not surprising that the many different jet-like 
phenomena that are observed in the solar atmosphere 
have been, and remain to some extent, such a puzzle. 
For example, in the same simulations, with similar mag- 
netic field geometry and driver, we find at different times 
jets that show significant chromosphcric signal without 
transition region or coronal counterparts, jets that show 
only blueshifts in the chromosphere, jets that only show 
redshifts in the chromosphere, jets that do not show chro- 
mosphcric signatures but are dominated by bidirectional 
flows in the transition region, and jets in the transition 
region lines that do not show any coronal counterpart. 
This is a natural consequence of the complex mix of mag- 
netic field geometry, the history of the plasma motions, 
the narrow height /temperature range in which the ob- 
servables are formed, and the varying mix of heating and 
acceleration in reconnection events. As a result, estab- 
lishing the relationship between various types of events 
in the solar atmosphere is fraught with difficulties that 
can only be resolved by statistical comparisons of spatio- 
temporal data of high quality with advanced radiative 
MHD models. 



This work was supported by the Research Council 
of Norway through grants 159483/V30, 170926, and 
170935, as well as a grant of computing time from the 
Program for Supercomputing. B. D. P. was supported 
by NASA grants NNG06GG79G, NNX08AL22G, and 
NNM07AA01C (Hinode). L. H. thanks Oystein Langan- 
gen and Mats Carlsson for useful discussions. 



Doyle, J. G., Popescu, M. D., & Taroyan, Y. 2006, A&A, 446, 327 
Erdelyi, R., De Pontieu, B., & Sarro, L. M. 1999, in ESA Special 

Publication, Vol. 448, Magnetic Fields and Solar Processes, ed. 

A. Wilson, 1345 
Grossmann-Doerth, U. & Schmidt, W. 1992, A&A, 264, 236 
Hansteen, V. H. 2005, in lAU Symposium 223, Multi- Wavelength 

Investigations of Solar Activity, ed. A. V. Stepanov, E. E. 

Benevolenskaya, & A. G. Kosovichev, 385 
Hansteen, V. H., De Pontieu, B., Rouppe van der Voort, L., van 

Noort, M., & Carlsson, M. 2006, ApJ, 647, L73 
Heggland, L., De Pontieu, B., & Hansteen, V. H. 2007, ApJ, 666, 

1277 

Innes, D. E. 2001, A&A, 378, 1067 



13 



Innes, D. E., Inhester, B., Axford, W. I., & Wilhelm, K. 1997, 

Nature, 386, 811 
Innes, D. E. & Toth, G. 1999, Sol. Phys., 185, 127 
Karpen, J. T., Antiochos, S. K., & DeVore, C. R. 1995, ApJ, 450, 

422 

Katsukawa, Y., Berger, T. E., Ichimoto, K., Lites, B. W., Nagata, 
S., Shimizu, T., Shine, R. A., Suematsu, Y., Tarbell, T. D., Title, 
A. M., & Tsuneta, S. 2007, Science, 318, 1594 

Madjarska, M. S. & Doyle, J. G. 2003, A&A, 403, 731 

Ning, Z., Innes, D. E., & Solanki, S. K. 2004, A&A, 419, 1141 

Osterbrock, D. E. 1961, ApJ, 134, 347 

Roussev, I., Doyle, J. G., Galsgaard, K., & Erdelyi, R. 2001a, A&A, 
380, 719 

Roussev, I., Galsgaard, K., Erdelyi, R., & Doyle, J. G. 2001b, A&A, 

370, 298 
— . 2001c, A&A, 375, 228 

Sarro, L. M., Erdelyi, R., Doyle, J. G., & Perez, M. E. 1999, A&A, 
351, 721 

Shibata, K., Nakamura, T., Matsumoto, T., Otsuji, K., Okamoto, 
T. J., Nishizuka, N., Kawatc, T., Watanabc, H., Nagata, S., 
UcNo, S., Kitai, R., Nozawa, S., Tsuneta, S., Suematsu, Y., 
Ichimoto, K., Shimizu, T., Katsukawa, Y., Tarbell, T. D., Berger, 
T. E., Lites, B. W., Shine, R. A., & Title, A. M. 2007, Science, 
318, 1591 



Shibata, K., Nishikawa, T., Kitai, R., & Suematsu, Y. 1982, 

Sol. Phys., 77, 121 
Sterling, A. C. 2000, Sol. Phys., 196, 79 

Sterling, A. C, Shibata, K., & Mariska, ,1. T. 1993, ApJ, 407, 778 
Suematsu, Y., Wang, H., & Zirin, H. 1995, ApJ, 450, 411 
Takeuchi, A. & Shibata, K. 2001a, ApJ, 546, L73 
— . 2001b, Earth, Planets and Space, 53, 605 

Tsiropoula, G., Ahssandrakis, C. E., & Schmieder, B. 1994, A&A, 
290, 285 

Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 
Wilhelm, K. 2000, A&A, 360, 351 



